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摘 要 : 太阳 系 外 类 地 行星 下 接 成 像 探 测 极 具 挑 战 性 。 在 可 见 光 短 波段 , 行星 信号 主要 
来 自 恒星 经 由 行星 大 气 的 反射 光 ,两 者 对 比 度 在 10” 量 级 ,导致 微弱 的 系 外 行星 光 被 淹没 。 
因此 ,有效 抑制 来 自 恒星 的 强 光 是 实现 类 地 行星 探测 的 关键 。 本 文 提出 了 基于 光 瞳 振幅 调制 
和 相位 校正 的 技术 方案 , 以 实现 全 区 域 高 对 比 度 成 像 。 针 对 上 述 方案 ,进行 了 数值 模拟 分 析 
来 说 明 其 可 行 性 及 潜在 性 能 。 数 值 模拟 基于 点 扩散 函数 的 能 量 分 布 来 评价 系统 的 目标 成 像 对 
HKE ,并 采用 随机 并 行 梯度 下 降 和 返 代 算法 进行 优化 ， 最 终 成 像 对 比 度 在 360 度 全 区 域内 获得 


达到 102 量 级 。 研 究 表 明 该 方案 有 望 用 于 未 来 空间 类 地 行星 直接 成 像 探 测 。 
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0 引 言 


太阳 系 外 行星 探测 是 当今 天 文学 领域 的 研究 热点 。 截 止 2019 年 5 月 ， 已 经 探测 发 现 4000 
多 颗 太 阳 系 外 行星 ”。 系 外 行星 探测 包括 间接 探测 和 直接 成 像 方法 。 间 接 方法 主要 包含 视 向 
速度 法 、 凌 星 法 天 体 测量 以 及 微 引 力 透镜 法 ， 是 通过 观测 恒星 物理 信息 反 演 ， 从 而 获得 相关 
的 行星 信息 ， 还 难以 直接 获得 行星 的 有 效 温度 、 大 气 等 重要 物理 信息 。 对 系 外 行星 进行 直接 
成 像 探 测 ， 能 够 测量 行星 的 质量 、 轨 道 ， 并 通过 进一步 光谱 分 析 能 够 研究 行星 大 气 组 成 、 表 
面 重力 和 有 效 温 度 等 重要 物理 信息 ， 是 未 来 确认 系 外 生命 特征 信号 的 关键 技术 ”。 

恒星 经 过 望远镜 孔径 会 产生 衍射 , 且 恒 星光 强度 远 高 于 行星 光 ， 导 致 微弱 的 行星 信号 被 
淹没 。 由 于 可 见 光 短波 段 存在 大 量 的 生命 特征 信号 ， 且 在 获取 同等 空间 分 辩 率 情况 下 (正比 
于 MD，》 为 工作 波长 ，D 为 望远镜 口径 ) ， 对 望远镜 口径 要 求 降低 。 因 此 ， 可 见 光 波段 是 
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来 成 像 搜 寻 类 地 行星 的 重要 工作 波段 0。 然 而 ， 在 短波 探测 围绕 类 太阳 光谱 型 恒星 宜 居 带 
内 的 类 地 行星 ， 成 像 对 比 度 需 要 达到 10-™P1。 高 对 比 度 星 网 仪 能 够 对 恒星 光 进行 有 效 抑制 以 
获得 超 高 对 比 度 成 像 区 域 。 近 年 来 , 国内 外 研究 学 者 们 提出 了 不 同 技术 原理 的 星 园 仪 , 理论 
ERGOT EGRET DAG BLOT, 但是， 星 罗 仪 系统 中 的 非 理想 光学 元 件 产 生 的 波 前 畸变 会 
带 来 散 斑 噪 声 ， 制 约 实际 的 成 像 对 比 度 。 在 没有 波 前 校正 的 情况 下 ， 成 像 对 比 度 只 有 10 E 
右 。 因 此 ， 有 效 较 正 波 前 畸变 以 抑制 消除 散 斑 噪声 是 实现 类 地 行星 超 高 对 比 度 成 像 的 关键 。 
Malbert 首 先 提 出 了 一 种 利用 可 变形 镜 (Deformable Mirror， 简 称 DM) 来 产生 暗 区 的 优化 控 
制 算 法 中， 该 方法 可 在 局 部 区 域 上 实现 高 达 10 的 成 像 对 比 度 。 

近 几 年 来 ， 得 益 于 可 变形 镜 技术 与 器 件 的 发 展 ， 星 冕 仪 实验 成 像 对 比 度 已 经 提高 至 10” 
量 级 。 然 而 ， 目 前 使 用 的 波 前 校正 方案 产生 的 高 对 比 度 成 像 区 域 较 小 ,无 法 实现 全 区 域 高 对 
比 度 成 像 ， 影 响 探测 探测 效率 VOM, 

针对 上 述 问题 ， 本 文 提出 一 种 360 度 全 区 域 、 超 高 对 比 度 成 像 星 冕 仪 技术 方案 ， 并 给 出 
初步 的 数值 模拟 结果 。 在 数值 模拟 过 程 中 , 采用 随机 并 行 梯度 算法 (Stochastic parallel gradient 
descent algorithm， 简 称 SPGD ) 进行 迭 代 优 化 ， 以 系统 PSF 能 量 反 馈 控 制 光 瞳 能 量 和 波 前 ， 
使 得 系统 成 像 对 比 度 达 到 最 优化 5 六。 整个 过 程 分 成 两 步 ， 首 先 通 过 优化 有 限 带 环 状 透 过 率 
渐变 滤 光 片 进行 振幅 调制 ， 在 4MD 处 获得 10 "的 对 比 度 ， 再 使 用 32x32 单 元 的 DM 进 行 波 前 校 
正 ， 最 终 在 成 像 焦 面 4-14MD 区 域内 获得 10” 的 高 对 比 度 成 像 暗 区 。 为 了 进一步 扩大 成 像 区 
ER, 还 提出 一 种 基于 振幅 调制 和 相位 校正 同时 优化 的 数值 模拟 方案 , 最 终 在 成 像 焦 面 4-20MD 
区 域内 获得 10" 的 高 对 比 度 成 像 暗 区 。 


1 优化 算法 原理 


恒星 光 统 通过 光学 系统 在 焦 面 处 的 点 扩散 函数 如 下 : 
Insp (x,y) = |F[P(u, v)e'?™”) |? (1) 
其 中 ，F 为 相关 函数 的 傅 里 叶 变 换 ，PGw 7) 为 光学 系统 入 瞳 处 的 振幅 ， 本 文 将 对 圆 形 孔径 进 
ITE; plu, 切 为 入 瞳 处 的 相位 。 星 冕 仪 系统 通过 在 瞳 面 处 优化 P(u 刀 和 9(wZ) 这 两 个 变 
量 来 获得 高 对 比 对 上 暗 区 ， 使 用 的 评价 函数 如 下 : 


= È lout y) 
J5 È lin (x,y) @) 
其 中 加 (cy 为 中 心 区 域 的 能 量 强度 ，1ouz(2 YVA RAR ER RD I HE, AVE DY eR 


于 找到 成 像 区 域 能 量 相对 于 中 心 区 域 能 量 的 最 小 值 。 本 文 使 用 的 高 对 比 成 像 区 域 范围 如 图 1 
所 示 ， 分 别 在 全 区 域 和 一 半 区 域内 进行 优化 。 
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图 1 定义 优化 区 域 示意 图 
Fig.1 Definition of optimization region 
成 像 区 域 的 对 比 度 定义 如 下 : 
I(x,y) 


max 


Contrast = 


(3) 
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其 中 ax 为 焦 面 处 的 能 量 最 大 值 ，1(x, y) 为 高 对 比 度 成 像 区 域内 某 一 点 的 能 量 。 根 据 圆 孔 点 
扩散 函数 的 分 布 特点 ， 其 能 量 主要 集中 在 中 心 艾 里 斑 内 ， 其 大 小 为 1.22MD，1nmax 为 艾 里 斑 
中 心 能 量 值 。 又 因为 圆 孔 点 扩散 函数 是 轴 对 称 的 , 本 文 将 选择 对 角 线 方向 来 评估 成 像 区 域 的 
对 比 度 。 
目标 函数 最 优化 有 诸多 不 同 的 优化 算法 ， 如 遗传 算法 由 、SPGD 算法 等 。 本 文 将 使 用 
SPGD 算法 对 振幅 和 相位 进行 迭代 优化 。 该 算法 通过 同时 施加 任意 小 的 扰动 到 所 有 的 变量 上 ， 
然后 评估 评价 函数 的 梯度 变化 ， 进 行 迭 代 优化 由。 变量 更 新 方程 如 下 : 

uktt = yk — yd] 8u" (4) 
其 中 是 迭代 次 数 ; = Uy, Uz a un 是 优化 变量 ; n 是 优化 的 参数 个 数 ; Y 表 时 增益 系数 〈 目 
标 优化 函数 最 小 时 取 正 ， 反 之 取 负 ) ; Su 为 服从 伯 努 利 分 布 的 任意 扰动 ，6J 是 评价 函数 的 
变化 量 : 


ôJ = J(utdéu) — J(u) = J(u, + buy,..., Un + un) — J us... Un) (5) 
为 了 提高 计算 精度 采用 双边 扰动 的 方法 ， 此 时 评价 函数 的 变化 量 为 : 
ôu ôu 
6) =I, -J-=J (u+) -J(u -F). (6) 
为 了 加 快 算法 的 收敛 速度 ,增益 系数 y 将 适应 于 评价 函数 变化 ， 此 时 第 £ 步 的 增益 系数 
7 为 : 


ere = yr JE (7) 


2 数值 模拟 


2. 1 振幅 调制 与 波 前 校正 分 步 优 化 数值 模拟 

星 冕 仪 结构 如 图 2 所 示 。 人 恒星 光 通 过 望远镜 焦点 (Telescope Focus) 并 由 准 直 透镜 Lı 

进行 准 直 ， 平 行 光 通过 孔径 光 阑 (Pupil Stop) 后 作用 在 有 限 带 环 状 透 过 率 渐变 滤 光 片 
(Stepped-transmission Filter) 上 进行 振幅 调制 ， 再 通过 成 像 镜 L 和 准 直 镜 Ls 在 暗 面 处 使 用 

DM 进行 波 前 校正 ， 最 后 通过 成 像 镜 Ly 进行 成 像 ， 在 成 像 焦 面 处 通过 探测 器 获 的 高 对 比 度 
HX 。 


Pupil Stop 


Telescope Focus Lens:L, Lens:L, DM 


Stepped-transmission filter 


Lens:L, 


Detector 


图 2 Ay RBS TA It A aS A 
Fig.2 Schematic diagram of the coronagraph with circular stepped-transmission filter 
使 用 SPGD 优化 算法 对 有 限 带 环 状 透 过 率 渐变 滤 光 片 进 行 优化 设计 。 滤 光 片 上 每 一 带 
的 透 过 率 作 为 一 个 独立 变量 ， 多 个 变量 进行 同时 和 迭代。 经 过 多 次 欠 代 优化 后 ,在 4-14MD 的 
区 域内 获得 107 的 对 比 度 。 图 3 为 优化 设计 的 有 限 带 环 状 透 过 率 渐变 滤 光 片 以 及 每 一 环 带 上 
对 应 的 透 过 率 。 图 4 为 经 过 有 限 带 环 状 透 过 率 渐变 滤 光 片 调制 后 在 焦 面 处 获得 PSF 图 像 。 
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Radius 


图 3 有 限 带 环 状 透 过 率 渐变 滤 光 片 示意 图 〈 左 ) 滤 光 片 沿 半径 方向 上 的 透 过 率 分 布 曲线 A) 


Fig.3 The transmittance distribution diagram of the stepped-transmission filter (Left) and the transmittance curve along the radius 


direction 
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图 4 使 用 有 限 带 环 状 透 过 率 渐变 滤 光 片 调制 后 ， 在 焦 面 获得 的 PSF 图 像 


Fig.4 The PSF image on focal plane modulated by stepped-transmission filter 


由 于 奈 硅 斯 特 采 样 定理 的 限制 "， 使 用 NxN 单元 的 DM 进行 波 前 校正 时 ， 其 可 调控 理 


论 效 范围 在 N/2 XAD 之 内 。 理 论 上 可 获得 的 成 像 对 比 度 可 如 下 经 验 公 式 表示 09: 


Contrast = 7 (=) 


2n0\2 
NA 


(8) 


其 中 co 为 镜面 RMS 值 ，N 为 DM 的 调 探 单元。 根据 式 〈8) 可 知 ， 可 以 通过 增加 DM 的 单元 


数 或 者 减 小 优化 范围 进一步 提高 对 比 度 。 


本 方案 模拟 了 一 个 32x32 的 可 变形 镜 ， H 


LE 论 优化 范围 为 16MD。 模拟 时 选用 理论 优化 范 


围 的 一 半 区 域 进 行 波 前 校正 。DM 中 每 一 块 子 镜 作为 一 个 独立 的 优化 变量 ， 经 过 多 次 迭代 优 


化 ， 在 4-8MD 的 全 区 域内 获得 3.5X102 的 对 比 度 ， 系 统 通过 率 32.04%。 图 5 为 DM 产生 
的 相位 分 布 图 以 及 经 过 DM 校正 之 后 在 焦 面 处 获得 PSF 图 像 〈 优 化 范围 4-8XD 全 区 域 )。 


图 6 为 4-8 和 ND 优化 区 域内 的 对 比 度 曲 线 ， 


中 黑色 实 线 为 未 调制 的 对 比 度 曲 线 ， 蓝 色 实 线 


~ 


为 振幅 调制 之 后 的 对 比 度 曲 线 , 红色 实 线 为 在 振幅 调制 的 基础 上 进行 波 前 校正 后 的 对 比 度 曲 


线 。 


图 5 DM 相位 分 布 图 (优化 范围 4-8MD)〉( 左 ) 使 


H DM 波 前 校 了 


1 
107 
o 10* 
10* 
10* 
10 


E 后 ， 在 焦 面 获得 的 PSF 图 像 〈 右 ) 


Fig.5 DM phase distribution diagram (4-8 MD optimization range) (Left) the PSF image on focal plane with DM wave-front correction 


(Left) 
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图 64-8 A/D 优化 


4 6 
Angular Separation (4/D) 


区 域内 对 比 度 曲线 


Fig.6 Contrast curve in the 4-8 MD optimization region 


在 4-14XMD 一 半 区 域 进行 波 前 校正 ， 经 过 
了 1.6X10" 的 对 比 度 ， 系 统 通过 率 32.04%. 


后 在 焦 面 处 获得 的 
对 比 度 曲 线 ， 其 中 


ee a 


VB UE ia» Æ 4-14/D 一 半 区 域内 获得 


PSF 图 像 〈 优 化 范围 4-14XD 一 半 区 域 )。 


图 7 DM 相位 分 布 图 〈 优 化 范围 4-14MD) CE) 使 


图 7 为 DM 产生 的 相位 分 布 以 经 DM 校正 之 


图 8 为 4-14MD 优化 区 域内 的 


eH EAT BL A EE 


DM 波 前 校正 后 ， 在 焦 面 获得 的 PSF 图 像 〈 右 ) 


Fig.7 DM phase distribution diagram (4-14 MD optimization range) (Left) the PSF image on focal plane with DM wave-front correction 


本 节 讨 论 进行 


在 两 个 调制 器 件 上 进行 同时 振幅 


(Right) 


Contrast 


Angular 


图 8 4-14XVD 优化 


Separation (MD) 


区 域内 对 比 度 曲线 


Fig.8 Contrast curve in the 4-14 MD optimization region 


2.2 振幅 调制 与 波 前 校正 同时 优化 数值 模拟 


同时 振幅 调制 和 波 前 相位 校正 。 恒 星光 通过 准 直 透 镜 Li 进行 准 直 后 作 月 


调制 和 波 前 校正 ， 然 后 通过 成 像 镜 Lo 进行 成 像 ， 在 成 像 焦 


ay 


pig 


面 上 通过 探测 器 获取 高 对 比 成 像 图 像 ， 其 结构 具体 如 图 9 所 示 。 
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图 9 AEA Se MR AE A TA IN OE SA AS ek 
Fig.9 Schematic diagram of the coronagraph with simultaneously optimizing amplitude and phase 
振幅 和 相位 在 均 为 环 带 分 布 , 每 一 坏 带 上 的 振幅 和 相位 为 一 个 独立 变量 , ZL TER 
优化 后 ,在 4-20MD 全 区 域内 得 到 对 比 度 为 1X10™" 的 成 像 暗 区 ,整个 系统 的 透 过 率 为 10.5%。 
10 为 产生 的 振幅 板 和 在 沿 半径 方向 上 的 振幅 分 布 曲线 。 图 11 为 产生 的 相位 板 以 及 经 振幅 
相位 优化 后 在 焦 面 处 获得 的 PSF 图 像 。 图 12 为 4-20XWD 优化 区 域内 的 对 比 度 曲 线 ， 其 中 红 
色 实 线 为 同时 振幅 调制 和 波 前 校正 后 的 对 比 度 aa 


pe os 
os os 
or o 
06 Eos 1 
2 | 
05 gos 
H 
oa 
04 = 
os 
o3 
02 1 
02 
oo 
ol 
N ll eae ) 
o so 190 150 
o 


图 10 振幅 分 布 图 〈 左 ) 治 半径 方向 上 的 透 过 率 分 布 曲线 〈 右 ) 


Fig.10 Amplitude distribution diagram (Left) and the transmittance curve along radius direction (Right) 


1 
2 
15 
1w? 
1 
os 
10* 
o 
25 
10% 
4 
A We 
2 
25 10 


图 11 相位 分 布 图 〈 左 ) 优化 后 在 焦 面 获得 的 PSF E E) 
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Fig.11 Phase distribution diagram (Left) the PSF image on focal plane after optimization (Right) 
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图 12 4-20 WD 优化 区 域内 对 比 度 曲线 


Fig.12 Contrast curve in the 4-20 MD optimization region 


3. 结论 


本 文 提出 的 两 种 优化 方案 在 全 区 域内 均 获 得 了 107 的 对 比 度 ， 为 今后 对 太阳 系 外 类 地 
J 星 直 接 成 像 提 供 了 可 能 。 表 1 为 三 个 数值 模拟 的 具体 结 
表 1 振幅 相位 分 步 优 化 和 振幅 相位 同时 优化 数值 模拟 结 


Tab.2 Results of phase optimization with setted amplitude and simultaneously optimization 


下 


amplitude and phase simulation 


Simulation 1 2 3 
IWA (A/D) 4 4 4 
OWA (A/D) 8 14 20 

Region Whole Half Whole 
Contast 3.5X10" 1.6X10™ 1x10" 
Throught (%) 34.04 34.04 10.5 


同时 优化 振幅 与 相位 可 以 在 较 大 区 域内 获得 更 高 的 对 比 度 , 该 方案 可 进行 大 区 域 全 视 场 
太阳 系 外 行星 探测 。 该 方案 的 透 过 率 较 低 ， 导 致 实际 观测 时 需要 较 长 的 曝光 时 间 ， 一 定 程 度 
上 影响 观测 效率 。 此 外 ， 同 时 优化 振幅 与 相位 所 需要 的 控制 单元 较 多 ， 目 前 DM 的 可 控 单 
元 数 无 法 满足 此 方案 的 需求 ， 实 际 实验 中 可 用 空间 光 调制 器 〈Spatial Light Modulator， 人 简称 
SLM) 进行 优化 调制 。 使 用 SLM 进行 振幅 调制 相 较 有 限 带 环 状 透 过 率 渐变 滤 光 片 ， 更 加 灵 
活 方便 ， 可 以 根据 暗 面 的 形状 进行 实时 调制 。 由 于 SLM 的 特性 ， 需 在 系统 中 加 入 两 块 的 线 
偏振 片 ， 导致 整个 系统 的 透 过 率 降低 。 下 一 步 将 在 目前 成 像 对 比 度 基 础 上 进行 透 过 率 优化 设 
计 ， 进 一 步 提 高 系统 透 过 率 ， 进 行 相关 的 测试 实验 ， 并 展开 宽 波 段 高 对 比 度 成 像 的 研究 。 
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Abstract: The direct imaging on Earth-like exoplanets faces great challenges. The faint 
planet light is the reflective light from the atmosphere of the planet. Thus, the intensity contrast 
between the planet and its host star is 107° in the visible wavelength, causing the planet light 
contaminated by the stronger stellar light. Therefore, it is critical to suppress the strong light from 
the host star in order to detect the Earth-like exoplanets. In this paper, we propose the apodized 
pupil amplitude modulation and phase correction technique to achieve an extreme high contrast. 


We then provide a numerical simulation to demonstrate the feasibility and potential performance. 


An iterative Stochastic Parallel Gradient Descent (SPGD) algorithm is applied in the numerical 
simulation to optimize the target imaging contrast by evaluating the point spread function. Finally, 
it has demonstrated a contrast on the order of 10™° in a whole 360-degree working area. The 
research demonstrates that the proposed approaches are promoting for the future Earth-like direct 
imaging detection in the space. 

Key words: Exoplanet; Terrestrial planet; Direct imaging; High contrast; Numerical 


simulation; Optimization algorithm 
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